Pairing-based cooling of Fermi gases 
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We propose a pairing-based metliod for cooling an atomic Fermi gas. A three component (labels 
1, 2, 3) mixture of Fermions is considered where the components 1 and 2 interact and, for instance, 
form pairs whereas the component 3 is in the normal state. For cooling, the components 2 and 3 
are coupled by an electromagnetic field. Since the quasiparticle distributions in the paired and in 
the normal states are different, the coupling leads to cooling of the normal state even when initially 
Tpaired > Tnormai- The cooliug efBcieucy is given by the pairing energy and by the linewidth of 
the coupling field. No superfluidity is required. The method has a conceptual analogy to cooling 
based on superconductor - normal metal (SN) tunneling junctions. Main differences arise from 
the exact momentum conservation in the case of the field-matter coupling vs. non-conservation of 
momentum in the solid state tunneling process. Moreover, the role of processes that relax the energy 
conservation requirement in the tunneling, e.g. thermal fluctuations of an external reservoir, is now 
played by the linewidth of the field. The proposed method should be experimentally feasible due to 
its close connection to RF-spectroscopy of ultracold gases which is already in use. 

PACS numbers: 03.75.Ss, 32.80.-t, 73.50.Lw 



Cooling is often a prerequisite for observing correlated 
quantum states of matter, as has been demonstrated, 
e.g., by the history of superconductivity, superfluidity 
and the quantum Hall effect. In the case of ultracold 
atoms, laser cooling and evaporative cooling have enabled 
the observation of Bose-Einstein condensates P, 0] , Mott 
insulator states ISl and condensates of strongly interact- 
ing Fermi gases [1 i, i, 0, H, SHE!. The research 
is rapidly expanding to new directions such as search for 
novel quantum states in optical lattices, p-wave pairing, 
and studies of multi-component mixtures of Fermions. To 
reach the required temperatures also in these new sys- 
tems, cooling that goes beyond the presently available 
methods may become necessary. We propose a method 
for cooling fermionic atoms in a normal state by an inter- 
nal tunneling contact to a paired state. The coexisting 
paired and normal states can be realized by trapping a 
three-component Fermi gas and tuning the various inter- 
component couplings by Feshbach resonances. The effec- 
tive tunneling contact is realized by coupling the internal 
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FIG. 1: (Color online) The normal state can be cooled by 
filling thermal holes (a) or removing hot particles (b). Net 
cooling is achieved because the thermal quasiparticle distri- 
butions are different in the (N) normal and (S) paired (en- 
ergy (pseudo)gap A) states. The picture is schematic: there 
is no spatial tunneling barrier, instead, the normal and paired 
states are coupled by an electromagnetic field. 



states of the atoms by electromagnetic fields. The cool- 
ing method is closely related to RF-spectroscopy of the 
pairing gap in ultracold Fermi gases l3, 12 1 and has an 
analogy to on-chip cooling methods of solid state struc- 
tures via superconductor - normal metal (SN) tunneling 
junctions [l3|. 

The basic principle of the proposed cooling is illus- 
trated in Fig. 1. A three-component Fermi gas is con- 
fined to a certain spatial volume, e.g. by an (optical) trap 
or lattice. Components 1 and 2 interact, for instance via 
a Feshbach resonance, whereas the component 3 does not 
interact with the other components. Components 2 and 
3 are coupled by an electromagnetic field, with a certain 
(effective) detuning 5 = lo field — ^23 and Rabi frequency 
n. All atoms coexist in the same spatial volume. Con- 
ceptually, however, one can imagine a paired state (1 and 
2) - normal state (3) interface where the electromagnetic 
field creates a tunneling coupling, as is displayed in Figs, 
la and lb to illustrate different processes that contribute 
to the cooling. In Fig. la, the detuning is chosen to be 
positive, S > 0, which means that the field gives to the 
system excess energy which is sufficient for breaking a 
pair: one component of the pair (1) becomes an excita- 
tion in the paired gas whereas the other one (2) is trans- 
ferred to the state 3 and fills a thermal hole in the normal 
gas. The inverse process, which also occurs for the same 
detuning, is that an atom from the normal state joins 
a thermal quasiparticle in the paired state and forms a 
pair. However, due to the pairing (pseudo)gap, there are 
fewer quasiparticles in the paired state than in the nor- 
mal state, see [ij] and discussion after Eq.Q, and this 
causes an imbalance between the two processes and thus 
a net cooling effect. In Fig. lb, a thermal quasiparticle 
in the paired state is transferred to the normal state. To 
enable this, the field has to take the excess energy of the 



2 



quasiparticle, therefore the process occurs for 5 < (in 
contrast to Fig. la where (5 > 0). Now the inverse pro- 
cess is the transfer of a hot quasiparticle away from the 
normal state. Again, the imbalance in the quasiparticle 
distributions leads to the desired coohng of the normal 
state. Note that the process corresponds to heat trans- 
fer between the paired and normal components, and no 
entropy is removed from the total system as will be dis- 
cussed in detail below. In the following, we demonstrate 
the efhciency of the method, determine by numerical and 
analytical calculations the optimal parameters and lim- 
itations for cooling, and discuss the connection and dif- 
ferences to SN tunneling junction coolers. 

We consider a three-component gas, i.e. fermionic 
atoms in three different internal states |2) and |3). 
Atoms in two of these states, |1) and |2), interact rea- 
sonably strongly. Components 1 and 2 can be different 
internal (e.g. hyperfine) states of an atom, or different 
atomic species (such as ^Li and ^°K). Components 2 and 
3 must be different internal states of the same atom since 
they have to be coupled by an electromagnetic field (for 
instance RF-field, or laser fields in a Raman configura- 
tion). Even if only pairing, not superfluidity, is necessary 
for the cooling, we assume here for convenience that the 
atoms in states |1) and |2) form a superfluid described 
by the standard BCS-theory with the BCS-Hamiltonian 
Hs- Atoms in the third state |3) are in the normal state, 
with the corresponding Hamiltonian Hn. Chemical po- 
tentials /Z5 and will be added to Hs and Hn, and 
incorporate also possible Hartree energies due to weak 
interactions between components |3) and |1) or |2), for 
definitions see [3, An electromagnetic field cou- 

ples the states |2) and |3). This atom-field coupling is 
described in the rotating wave approximation, and the 
total Hamiltonian becomes 



H = Hn + Hs + -{Ns - Nn) + Ht. 



(1) 



The atom-field coupling (tunneling) part of the Hamilto- 



nian is Ht — X]fc ^ '^k 



Vt* clj^Cks {h = 1). We treat Ht as a perturbation and 
calculate the rate of change of the number of atoms in the 
superfiuid state, (Ns) = -{Nn) = j-{[H,Ns]) = (N)^ 
the lowest order, i.e. using linear response theory [14| . 
From this we obtain the number of transferred atoms at 
time T as (A7V(r)) = J^{N{f))df. The delta-functions 
which enforce energy conservation in (N) are approxi- 
mated by (Cauchy-)Lorentz distributions of the width F, 
C(x) — — -i^T^-i . We assume that the intrinsic linewidth 
of the transition is small and the total linewidth is essen- 
tially given by the pulse Rabi frequency, or inverse pulse 
length, i.e. F ^ $7. (The results are not sensitive to the 
lincshape as shown in the following by comparing results 
for Gaussian and Lorentz distributions). We set F = 17, 
and fir = tt corresponding to a 7r-pulse which maximizes 



the transfer. This leads to 
{AN{t = ir/n)) = 27r^ 

k 

fi2 



{Ek-ii + sr + n^' 
(Ek + ^i-sr + n^ 



IKiv(ef)-«F.s(i?fe))+ (2) 
vl {nFMC)-nFA-Ek)). 



Here ^"^/^ - ^ 



fiN/s and E, = ^iCiy + A^. The 
number of thermal quasiparticles is given by the Fermi 
distribution \/{e^/^'^^ -f 1); in the paired state E = Ek 
with the minimum value A whereas in the normal state 
E — ^j^ and there is no gap, the distributions in the 
paired and normal states are therefore different. In the 
following, unless stated otherwise, we assume fiN = fJ'S, 
Tjv = Ts and set h = ep = fcs = 1- For clarity, we keep 
F and Q explicit although in the numerics we set F = Jl. 
The Eq.([2|) shows exphcitly the processes illustrated in 
Fig.l. The term proportional to Uk corresponds to the 
single quasiparticle transfer in Fig. lb, and the term pro- 
portional to Vk to the pair breaking/formation process in 
Fig. la. The differences in the quasiparticle distributions, 
given by np^is[{£^j^ ) and np^siEk), determine the direction 
of the net particle transfer. For pairing gap A = 0, there 
is no net particle transfer and no net cooling. 

We also derive the heat fiux out of the normal state (i.e. 
the cooling power) from (Qn) = {NHj^) in the same way 
as (N) above. Since Hn measures energy from the chem- 
ical potential hn, the energy of removed (inserted) hot 
quasiparticles enters with the same sign as inserted (re- 
moved) particles that fill (create) holes below the Fermi 
level. Therefore (Qn) directly tells about the amount of 
cooling. In the case of solid state SN junctions, the heat 
flux Q is calculated in the same way but the functional 
dependence on the relevant parameters differs due to the 
momentum conservation in our case. We obtain 

k ^ 

4 C {^1 -Ek-S) [np,Ni^^) - npMEk)] + (3) 
vl C (Cf + Ek-S) [np^Ni^l^) - nF,s{-Ek)] 

To determine the final temperature of the normal state 
gas after the cooling pulse, we assume thermalization 
into a Fermi distribution after the pulse (for instance 
due to weak p-wave interactions), or that the distribu- 
tion produced by the cooling pulse is close to a Fermi 
distribution. We then obtain the final temperature T 
and the chemical potential jl from (Nn) + (AN) = 
Efe«F(lfc) and (En) + {AE) = Y^k ^"^(6), where 
En — np is the Fermi function with 

T and jl, and {AE) is the energy added by transferring 
the (AA^) particles. 
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FIG. 2; (Color online) Temperature of the normal state af- 
ter the cooling pulse as a function of detuning, for different 
linewidths V = 0.01, 0.05, 0.16 and O.SOef. Before the pulse, 
the temperature was 0.2Tf and the pairing gap is 0.2eF. The 
linewidth 0.16 ± O.Olep gives the maximal cooling. 



FIG. 3: (Color online) Energy relative to chemical potential 
transferred out of the normal state (heat flux), (Qjv), as a 
function of the detuning S. The curves are numerical results 
for Gaussian and Lorentz distributions, and approximate ana- 
lytical result for a Gaussian distribution. The gap A — OAep, 
T = 0.1Tf and the linewidth F = O.OSef- 



be in the limits \d\ < A. The result becomes 



In Fig. 2 is shown the temperature after the cooling 
pulse, as a function of the detuning, for four different 
linewidths F. The optimal detuning 6 is somewhat be- 
low the pairing gap energy (in general we find S ~ A/2), 
and the optimal linewidth F('^ 17) is of the same order 
of magnitude. We have tested also temperatures much 
smaller than the gap: then there is nearly no transfer of 
particles below the gap energy due to the lack of ther- 
mal quasiparticles, moreover, the cooling effect is minute 
and only heating is caused for larger detunings. In other 
words, the pairing gap and the temperature should be 
of the same order of magnitude, A ~ T. In the limits 
T >> A and T << A the desired effect vanishes. 



Note that for larger linewidths, the cooling works also 
for zero detuning S. Since S plays the same role as volt- 
age applied over an SN contact, this corresponds to cool- 
ing in presence of no external voltage. Such a system 



has been discussed recently in an interesting way in 17 1 
where cooling by an SN junction is driven solely based on 
thermal noise from a resistor - the system is viewed as 
a Brownian refrigerator which conveys heat unidirection- 
ally in response to external noise, in analogy to Brow- 
nian motors and thermal ratchets or to the concept of 
Maxwell's demon. The effect of the thermal noise is to re- 
lax the energy conservation requirement in the tunneling 
process - in our case, this is done by the finite linewidth 
F of the coupling field. We have derived an analytical 
estimate for the heat flux (Qn) defined in Eq.(l3]), using 
the following approximations: the Fermi distributions are 
approximated by their Boltzmann-like tails, the Lorentz 
distribution is replaced by a Gaussian of the same width, 
and the temperature and the linewidth are assumed to 
be smaller than the energy gap, also the detuning should 
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Here T = Tj^. The result can be compared to the one in 



17 1 by substituting ^ a = 2TrEc and S E^. In 
17| the variance a for the voltage arises from the temper- 
ature of the environment Tr, and contains the charging 
energy Ec due to the junction and other capacitances. 
All the exponential factors are the same in our result 
and [l3|, however, the other terms differ, originating from 
the momentum conservation in our case. The analytical 
estimate (|H) is compared to numerical results in Fig. 3 
by plotting the numerically calculated heat fluxes for a 
Gaussian and a Lorentz distribution, and the approxi- 
mated result ([3]) for a Gaussian distribution. Fig. 3 shows 
that the analytical estimate is reasonably good around 
the optimal detuning for cooling, however, one should use 
it with care for (5 ~ 0. From Eq.(I3]) one can also estimate 
the optimal value for the linewidth, in two ways: one can 
either require a) the argument of the exponential to be 
zero, or b) those terms in the prefactor that are nega- 
tive to be minimized. This gives r = yJxT{A + 5) where 
X — 2 for case a) and x = 1 for case b). Notably, this 
is the same condition as obtained in (l7j for the optimal 
amount of "energy uncertainty" (variance a oc Tf?) in the 
cooling process. For the parameters of Fig. 2 the optimal 
linewidth given by these estimates is between Q.2ep (a) 
and QAAep (b) which is in excellent agreement with the 
numerically found optimal value 0.16 ± Q.Qlep. 
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Similarly, the heat flux out of the paired state becomes 

Qs _ STTil^VJ^ 



N 



N 



(2A- 



T, 



N 



(5) 



where one can approximate the last term to be zero when 
r^/TV << A. One can see from here that, for suitable 
values of < T5, it is possible to cool the paired state 
as well. However, if Tjv > Ts cooling of the paired state 
is not possible because the factor ^ + -^^r is always 

positive. This is in contrast to the case of where 
the superconductor can be cooled (and the normal metal 
heated) also for Tn = Ts in the case that the external re- 
sistor connected to the SN junction has a temperature Tji 
lower than Tn- This illuminates the difference between 
the external reservoir (resistor) and the electromagnetic 
field: although thermal fluctuations of the reservoir in 
ITj l and the linewidth of the field in our case have the 
same role in relaxing the energy conservation in tunnel- 
ing, the reservoir can also increase its entropy in the pro- 
cess (if Tr < T/v) whereas the electromagnetic field does 
not have a role in the trade of entropy. In our case the 
entropy decreased by cooling the normal state simply in- 
creases the entropy of the paired state, and vice versa. 
E.g. if normal state is cooled from 0.2 Tp to 0.1 Tp the 
paired state heats from 0.2 Tp to 0.29 Tp. To continue 
the cooling efficiently, the hot quasiparticles produced by 
the cooling pulse could be selectively removed away from 
the paired state, to a state that is not confined in the 
trap/lattice. 

The proposed cooling method is not sensitive to hav- 
ing the chemical potentials equal, 113 — fJ-N, nor to homo- 
geneity or completeness of the gap. This is demonstrated 
in Fig. 4 where we show the heat flux for a non-balanced 
case fiN 7^ Ms: cooling works even better, it only be- 
comes asymmetric with respect to the detuning S. In 
practice, the paired state may well be different from the 
homogeneous space BCS-state used here, for instance the 
(pseudo)gap may be non-complete (i.e. it may have some 
spectral weight at zero energy as well), x-dependent, k- 
dependent, etc. In Fig. 4 we show a linear combination of 
results for different gaps, simulating nodes and inhomo- 
geneities in the gap in the spirit of local density approx- 
imation. The cooling is not dramatically affected, it is 
simply reduced in proportion with the inhomogeneity of 
the gap. The scheme also works for molecules but is not 
optimal because the pairing energy in that case easily be- 
comes much larger than the temperature. In general, the 
cooling is based on coupling systems with different spec- 
tral densities which can be caused also by other effects 
than pairing, e.g. by Landau quantization [lit . 

The setting described here can be used for ther- 
mometry as well, as is also done in the solid state 
context [l3| . The number of transferred particles be- 



c 
3 



X 

3 



ra 



-0.4 




FIG. 4: (Color online) The heat flux {Qn} out of the normal 
state. The curves are for various chemical potential difi'er- 
ences A/i = /is — /iat between the paired and the normal state. 
The gap A = 0.3eF, T = O.ITf and the linewidth O.OSeF. A 
linear combination of heat fluxes with different energy gaps, 
Qn = 0.1Q]v(A = 0) + 0.3Qiv(A = 0.15) + 0.6Qjv(A = 0.3), 
is also shown. 



comes exponentially sensitive to temperature; calculat- 
ing it in the same way as above gives the analyti- 



cal estimate 
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In summary, we have proposed cooling of a normal 
state Fermi gas, based on paired state - normal state 
tunneling created by electromagnetic fields. We have de- 
termined the optimal parameters for cooling and shown 
that temperature drops of the order of the pairing gap 
can be achieved. The cooling of a normal Fermi gas may 
find several applications. A promising one is to produce a 
very cold single-component Fermi gas, and subsequently 
turn on p-wave interactions by Feshbach resonance tech- 
niques. Especially in optical lattices, studies of p-wave 
pairing are of great interest due to the richness of the 
phase diagram. Our work also opens up a new con- 
nection and an opportunity for synergy between ultra- 
cold gases and solid state nanostructures research. In 
general, cooling-techniques are related to manipulation 
of the quasiparticle distributions of a quantum system, 
which is essential in the study of non-equilibrium effects. 
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